QED corrections to the parity-nonconserving 6s-7s amplitude in 133 Cs 
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The complete gauge-invariant set of the one-loop QED corrections to the parity-nonconserving 6s-ls ampli- 
tude in 133 Cs is evaluated to all orders in aZ using a local version of the Dirac-Hartree-Fock potential. The 
calculations are peformed in both length and velocity gauges for the absorbed photon. The total binding QED 
correction is found to be -0.27(3)%, which differs from previous evaluations of this effect. The weak charge of 
133 Cs, derived using two most accurate values of the vector transition polarizability (3, is Qw = —72.57(46) 
for 13 = 26.957(51)a| and Q w = -73.09(54) for (3 = 27.15(ll)o|. The first value deviates by 1.1a from 
the prediction of the Standard Model, while the second one is in perfect agreement with it. 
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PACS numbers: 1 1.30.Er,31.30.Jv,32.80.Ys 

Investigations of parity noncoservation (PNC) effects in 
atomic systems play a prominent role in tests of the Stan- 
dard Model (SM) and impose constraints on physics beyond it 
1 1, 2]. The 6,s-7s PNC amplitude in 133 Cs |3] remains one of 
the most attractive subject for such investigations. The mea- 
surement of this amplitude to a 0.3% accuracy 0.S has stim- 
ulated a reanalysis of related theoretical contributions. First, 
it was found o IE 01 that the role of the Breit interaction 
had been underestimated in previous evaluations of this effect 
11(1 ll ill . Then, it was pointed out 11211 that the QED cor- 
rections may be comparable with the Breit corrections. The 
numerical evaluation of the vacuum-polarization (VP) correc- 
tion 1 13] led to a 0.4% increase of the 6s-7s PNC amplitude in 
133 Cs, which resulted in a 2.2cr deviation of the weak charge 
of 133 Cs from the SM prediction. This has triggered a great 
interest to calculations of the one-loop QED corrections to the 
PNC amplitude. 

While the VP contribution can easily be evaluated to a high 
accuracy within the Uehling approximation, the calculation of 
the self-energy (SE) contribution is a much more demanding 
problem (here and below we imply that the SE term embraces 
all one-loop vertex diagrams as well). To zeroth order in aZ, 
it was derived in Refs. fl4lO . This correction, whose rel- 
ative value equals to — a/(27r), is commonly included in the 
definition of the nuclear weak charge. The aZ-dependent part 
of the SE correction to the PNC matrix element between s 
and p states was evaluated in Refs. Ilq.ll7ll . These calcula- 
tions, which are exact to first order in aZ and partially include 
higher-order binding effects, yield the correction of -0.9(1)% 
OH] and -0.85% 1 17]. This restored the agreement with 
SM. 

Despite of the close agreement of the results obtained in 
Refs. 00, the status of the QED correction to PNC in 
133 Cs cannot be considered as resolved until a complete aZ- 
dependence calculation of the SE correction to the 6s-7s tran- 
sition amplitude is accomplished. The reasons for that are the 
following. First, in case of cesium (Z = 55) the parameter 
aZ as 0.4 is not small and, therefore, the higher-order cor- 
rections can be significant. Second, because the calculations 

EES El 

are performed for the PNC matrix element only, 



they do not include other SE diagrams which contribute to the 
6s-7s transition amplitude. For instance, these calculations do 
not account for diagrams in which the virtual photon embraces 
both the weak interaction and the absorbed photon. Our calcu- 
lations, however, show that the contributions of all diagrams 
are of the same order of magnitude (in both length and ve- 
locity gauges, see below), and the final result arises through a 
delicate cancellation of individual terms, none of which can be 
neglected. Third, strictly speaking, the PNC matrix element 
between the states of different energies is not gauge invariant. 
Despite the gauge-dependent part is suppressed by the small 
energy difference I17II . estimates of the uncertainty in the def- 
inition of the PNC diagrams may fail due to unphysical origin 
of the gauge-dependent terms. 

The first step towards a complete aZ-dependence calcula- 
tion was done in Ref. 1 19], where the SE correction to the 2s- 
2pi/2 PNC transition in H-like ions was evaluated. This tran- 
sition was chosen to deal with the simplified gauge-invariant 
amplitude. The results of that work agree with those of Refs. 
Ilall7lll 811 . However, as was stressed there, no claims can be 
made about the applicability of these results to the 6s-7,s PNC 
transition in neutral cesium. In this Letter we calculate the 
whole gauge-invariant set of the one-loop QED corrections to 
the 6s-7,s PNC transition amplitude in 133 Cs and compare the 
obtained result with the previous treatments. 

A systematic derivation of the QED corrections in a fully 
relativistic approach requires the use of perturbation theory 
starting with a one-electron approximation in an effective lo- 
cal potential V(r). In neutral atoms, it is natural to assume 
that V(r) includes not only the Coulomb field of the nucleus 
but also a part of the electron-electron interaction. The in- 
teraction of the electrons with the quantized electromagnetic 
field and the correlation effects are accounted for by perturba- 
tion theory. In this way we obtain quantum electrodynamics 
in the Furry picture. 

To derive formal expressions for the transition amplitude 
we employ the method developed in Ref. 12011 and described 
in detail in Ref. |21]. Since the wave length of the absorbed 
photon is much larger than the atomic size, one can use the 
dipole approximation. Within this approximation, calcula- 
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tions in the velocity gauge are performed using formulas given 
in Ref . 1 2 1 ] with the replacement exp (ik • x) — > 1 in the pho- 
ton wave function. The corresponding formulas in the length 
gauge are easily obtained by replacing a. with r in all vertices 
corresponding to the photon absorbtion and by multiplying the 
amplitude with the factor i(E b — E a ), where E a and E b are 
the total energies of the atom in the initial (6s) and final (7s) 
states, respectively. This simple rule can be derived using Eq. 
(205) of Ref. I2II1 and the equal-time commutation relations. 

To zeroth order, the 6s-7s PNC transition amplitude, which 
is usually employed in these calculations, is 
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Here a and b denote the 6s and 7s one-electron states, respec- 
tively, with the angular momentum projections m a = m b — 
1/2, d z — ez is the z projection of the dipole moment opera- 
tor (e < 0), Hw = — {G f / ^/8)Qw Pnuc{r)j5 is the nuclear 
spin-independent weak-interaction Hamiltonian Jll, Gf is the 
Fermi constant, 75 is the Dirac matrix, and p nuc is the weak- 
charge distribution. The one-loop SE corrections are defined 
by diagrams presented in Fig. [2 The derivation of the for- 
mulas for these diagrams is very similar to that for the QED 
corrections to the transition amplitude described in detail in 
Ref. 12 ill . As a result, the SE correction is given by the sum 
of the following terms: 
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Here the SE operator is defined by 
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I{uj) = e 2 a> i a v D liV {u)), = 7 V = (La), 
is the photon propagator defined as in Ref. |21], E'(-B) = 
dYj(E)/dE, and u = 1 — iO ensures the correct position 
of poles of the electron propagators with respect to the in- 
tegration contour. Taking into account the corresponding di- 
agrams with the mass counterterm results in the replacement 
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TABLE I: The binding energies of low-lying states in Cs, in a.u. 



State 


Loc. pot. 


DHF 


Exp. 


6S1/2 


-0.13079 


-0.12824 


-0.14310 


6P1/2 


-0.08696 


-0.08582 


-0.09217 


7S1/2 


-0.05621 


-0.05537 


-0.05865 


7Pl/2 


-0.04251 


-0.04209 


-0.04393 



T,(E) — > £(2?) — 7°i5m. The expressions for the VP correc- 
tions, which do not contain any insertions with the external 
photon line or the weak interaction attached to the electron 
loop, are obtained from Eqs. Q-Q by the replacement of 
the SE operator with the VP potential. The other VP correc- 
tions will not be considered here, since their contribution is 
negligible. 

The corresponding expressions in the velocity gauge are ob- 
tained by the replacement d z — ► —iea z / {E^ — E a ), where the 
energies E a and Eh include the QED corrections. In addition 
to the replacement d z — > —iea z /(eb — s a ) in Eqs. (fTl- (ll3> . it 
yields the contribution 

c™dd (&|S( £b )|&)-(q|E( £a )|q) 

which results from the expansion of the factor 1 / (Eb — E a ). 

Formulas ||2V(I15> contain ultraviolet and infrared diver- 
gences. To cancel the ultraviolet divergences, we expand 
contributions (0-0 mto zero-, one-, and many-potential 
terms and contributions ll8l- (ll 1> into zero- and many-potential 
terms. The ultraviolet divergencies are present only in the 
zero- and one-potential terms. They are removed analytically 
by calculating these terms in the momentum space (for details, 
we refer to Refs. j^H^Hil)- The many -potential terms are 
evaluated in configuration space. The infrared divergences, 
which occur in contributions (0-l|5j and (I12M 13>. are regu- 
larized by introducing a nonzero photon mass and cancelled 
analytically. 

Since the levels 6s, 6pi/2, 7s, and 7pi/2 are very close to 
each other, to get reliable results for the transition amplitude 
under consideration, one needs to use a local potential V(r) 
that reproduces energies and wave functions of these states to 
a sufficient accuracy. We construct such a potential by invert- 
ing the radial Dirac equation with the radial wave function 
obtained by solving the Dirac-Hartree-Fock (DHF) equation 
with the code of Ref . l25ll . Details of this procedure will be 
published elsewhere. In Table 1, we compare the energies ob- 
tained with the local potential V(r), that was derived using 
mainly the DHF wave function of the 6s state, with the DHF 
energies and with the experimental ones. 

Numerical evaluation of expressions (fH- (ll5> was per- 
formed by employing the dual-kinetic -balance finite basis set 
method 1 26] with basis functions construced from B-splines. 
The calculation of the zeroth-order contribution, with V(r) 
constructed as indicated above, yields £pnc =-1-002, in units 
ixlO _11 <5vK/(~^) a - u - This value should be compared with 
the corresponding DHF value, -0.742 |7], and with the value 



that includes the correlation effects, -0.908 OH. The results 
for the SE corrections are presented in Table 2. Since there is a 
significant cancellation between terms containing the infrared 
singularities, the terms corresponding to n = a in S'(e a ) 
and n = b in T,'(e a ) are subtracted from contributions (0- 
(0 and added to contributions J12i - (I13I . The total SE cor- 
rection <5-Ep°j!j c contains also the free term, — a/(27r)£'pNc^ 
mentioned above. Since this term is usually included into the 
weak charge Qw> one has to consider the binding SE correc- 
tion denned as SE^ = 5Ep& c + a/(2ir)E PNC - Accord- 
ing to Table 2, the binding SE correction amounts to -0.67%. 
To estimate the uncertainty of this value due to correlation 
effects, we have also performed the calculations with V(r) 
constructed employng the DHF wave function of the 7s state. 
While this leads to a 2% decrease of the transition amplitude, 
the relative shift of the SE correction is five times smaller. 
Since the correlation effects contribute to the transition am- 
plitude on the 20% level, we assume a 4% uncertainty for the 
total SE correction. Therefore, our value for the binding SE 
correction is -0.67(3)%. This value differs from the previous 
evaluations of the SE effect, -0.9(1)% J3 and -0.85% 0. 

We have also calculated the VP correction. Our value for 
the Uehling part amounts to 0.410%, which agrees well with 
the previous calculations of this effect. We have found that 
including the screening into the Uehling potential does not 
affect this value. As to the Wichmann-Kroll (WK) correction, 
our calculation employing approximate formulas for the WK 
potential from Ref. O yields -0.004% (cf. |9]). This leads 
to the 0.406% result for the total VP correction. Therefore, 
the total binding QED correction amounts to -0.27(3)%. 

To get the total 6s-7s PNC transition amplitude in 133 Cs, 
we combine the value that includes the correlation effects 
-0.908(1)± 0.5%, with the -0.61% Breit correc- 
tion the -0.27(3)% binding QED correction, the -0. 19(6)% 
neutron skin correction 1 28], the -0.08% correction due to the 
renormalization of Qw from the atomic momentum transfer 
q — 30 MeV down to q = 0, and the 0.04% contribution 
from the electron-electron weak interaction 1 17]. Using the 
experimental value for the Ep^c/P 01- where (3 is the vec- 
tor transition polarizabilty, we obtain for the weak charge of 
133 Cs: 

Qw = -72.57(29) cxp (36) th (16) 

for/3 = 26.957(51)a| HH and 

Q w = -73.09(39) oxp (37) th (17) 

for (3 = 27.15(ll)a 3 3 USUI. We conclude that the first 
value deviates from the SM prediction of -73.09(3) fiHl by 
1.1. a, while the second one is in perfect agreement with it. 

In summary, we have calculated the QED correction to the 
6s-7s PNC transition amplitude in 133 Cs and derived the weak 
charge using two most accurate values of the vector transi- 
tion polarizability. Further improvement of atomic tests of 
the Standard Model can be achieved, from theoretical side, by 
more accurate calculations of the electron-correlation effects 
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TABLE II: The SE corrections to the 6s - 7s PNC amplitude in 
133 Cs, in %. The results are presented in both the length (L) and the 
velocity (V) gauge. 



Contr. L-gauge V-gauge Contr. L-gauge V-gauge 

8E£ NC 4109 4X11 5E£ NC A04 -3.40 

&Ep N c 131 1.11 (S^pnc -4.61 -3.97 

5£p NC 0.34 0.40 SE i PNC 1.49 1.73 

(5£pnc -0-38 -0.32 SE^c -0.79 -1.03 

<5£p NC -1.29 -1.53 <S£ PNC 2.05 1.41 

<5£ PNC 3.89 3.25 SE^ d c 0.00 0.10 

5Ef, NC 1.33 1.57 SE^c -0.79 -0.79 

6E$S& -0.67 -0.67 



and, from experimental side, by more precise measurements 
of the PNC amplitude in cesium or other atomic systems. Par- 
ticularly interesting is the francium atom, where PNC effects 
are greatly enhanced by strong electric field of the nucleus. 
Precise measurements of the PNC amplitude in Fr are becom- 
ing feasible due to recent advances in producing, storing, and 
cooling of short-lived radioactive isotopes. 
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FIG. 1 : Feynman diagrams for the SE corrections to the PNC transition amplitude. The wavy line terminated with a triangle 
indicates the absorbed photon. The dashed line terminated with a cross indicates the electron-nucleus weak interaction. 



